Magnetic phases in the correlated Kondo-lattice model 
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We study magnetic ordering of an extended Kondo-lattice model including an additional on-site 
Coulomb interaction between the itinerant states. The model is solved in the dynamical mean-field 
theory using Wilson's numerical renormalization group approach as impurity solver. For a bipartite 
lattice we find at half filling the expected antiferromagnetic phase. Upon doping this phase is 
gradually suppressed and hints towards phase separation are observed. For large doping the model 
exhibits ferromagnetism, the appearance of which can at first sight be explained by Rudermann- 
Kittel-Kasuya-Yosida interaction. However, for large values of the Kondo coupling J significant 
differences to a simple Rudermann-Kittel-Kasuya-Yosida picture can be found. We furthermore 
observe signs of quantum critical points for antiferromagnetic Kondo coupling between the local 
spins and band states. 
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I. INTRODUCTION 

The Kondo-lattice model (KLM) is one of the most 
analyzed models in solid state theory due to its large 
variety of applications. Quite generally, the KLM de- 
scribes the interaction between localized spins and a band 
of conduction electrons. One particular class of sys- 
tems, where such a situation can be realized, are the 
transition metal oxides, which have received consider- 
able interest due to their rich phase diagram comprising 
different types of magnetic and orbital order and even 
supercondunctivityi This complexity stems from the in- 
terplay between the formation of narrow 3e?-bands lead- 
ing to a derealization of these states on the one hand and 
the local part of the Coulomb interaction between the 
3d-electrons tending to localize themi A particularly in- 
teresting example is Lai-^Ca^MnO^. In this cubic per- 
ovskite the five-fold degenerate 3d level is split by crystal 
field into three- fold degenerate t2 g , which have the lower 
energy, and two-fold degenerate e g states. These states 
have to be filled with 4 — x electrons, nominally yielding 
a metal even for x = 0. However, taking into account 
the local Coulomb interaction, three of these electrons 
will occupy the f^g-states forming an S = 3/2 high-spin 
state due to Hund's coupling, which interacts ferromag- 
netically with the electron occupying the e g states. Be- 
sides its complicated phase diagram with a large variety 
of paramagnetic and magnetically ordered metallic and 
insulating phases one finds a colossal magneto-resistance 
(CMR)4 The physics just described can be covered in 
great parts by the KLM, in this context usually called 
Double Exchange Model (DE)4£ 

Another class of materials which can be addressed 
by the KLM with ferromagnetic exchange interaction 
between conduction states and the localized spin de- 
grees of freedom are magnetic semiconductors or semi- 
metals in the series of the rare earth monopnictides and 
monochalcogenidesi£i£ Here the focus lies on the mag- 
netic respectively magneto-optic properties of for exam- 
ple EuS, EuO, GdN or thin films of such systems £ 

Last but not least the ferromagnetic KLM can be ap- 



plied to investigate the magnetic properties of diluted 
magnetic semiconductors such as Gai-^Mn^As. In these 
materials the III-V semiconductor GaAs becomes ferro- 
magnetic by introducing magnetic ionS ) 10 i n such as Mn. 
The local spin of the Mn-ions couples ferromagnetically 
to the states of the semiconductor, a setup which can be 
modeled by the KLM. Note, however, that in these ma- 
terials disorder effects are expected to play a vital role^ 

A completely different realization of the KLM starts 
from the periodic Anderson model (PAM) ; 13 i 14 i 15 a 
model for heavy fermion physics^ Heavy fermion 
physics manifests itself in a number of lanthanide and 
actanidc compounds, which have a very large effec- 
tive mass in common. The low-temperature physics of 
these compounds is determined by a partly filled f shell 
and hybridization induced spin flip scattering between 
the f and conduction electrons. By a Schrieffer-Wolff- 
Transformationii the periodic Anderson model can be 
mapped onto the KLM with antiferromagnetic exchange 
interaction. Note that in these materials one typically ex- 
pects a competition between the heavy fermion physics, 
driven by the Kondo effect due to the antiferromagnetic 
coupling, and the formation of antiferromagentism.— 

In most of the approaches one does not include 
Coulomb interaction among the conduction band elec- 
trons. Especially for manganitcs this is quite likely an 
insufficient approximation, because all the physics takes 
place in a correlated d-band, as explained above. There 
is no reason to ignore the local Coulomb correlations in 
the itinerant subshells, in particular because estimates 
of its magnitude typically lead to values of the order or 
even larger than the bandwidth of the d states at the 
Fermi level X Therefore, we want to address the influence 
of these local Coulomb correlations among the conduc- 
tion band on the magnetic properties of the KLM. 

Generally, the 5-fold or 14-fold degenerate d- respec- 
tively f-shells split in a crystalline environment into more 
or less well separated subshells. Here, we assume that 
the different crystal field levels are well separated on the 
scale of relevant low-energy structures - this is typically 
true in transition metal compounds, but less obvious in 
rare earth systems - and focus on a situation, where 
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we have one subshell well localized, hosting a spin S, 
while the remaining levels are split sufficiently to leave 
one relevant, spin-degenerate state close to the Fermi en- 
ergy. Our model thus consists of a conventional one-band 
Hubbard-Model ) 19 ' 20 i 21 where the itinerant states are in 
addition coupled to local spins S, i.e. 

H = i?Hub + -f^Spin ■ (1) 

Hnub describes the ordinary Hubbard Model, 



Hub 
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where cj CT (ci j(T ) creates (annihilates) an electron at lattice 
site Ri with spin a. While the inclusion of an arbitrary 
hopping tij presents no principle problem, we focus on 
the simplest case of next-neighbor-hopping parametrized 
by an amplitude t. Finally, hi^ a = cf a Ci, a denotes 
the density operator for a a— electron at site Ri and U 
parametrizes the local Coulomb interaction. 

The band electrons are in addition coupled to a local 
spin Si at each lattice site by an exchange interaction 
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J ^ ^ Si • Si 



where Si is the spin operator for the band states at site Ri. 
As discussed before, such a term can arise through Hund 
coupling, in which case J is ferromagnetic, or through a 
hybridization, leading to an antiferromagnetic J Note 
that both effects can appear simultaneously, thus par- 
tially compensating each other.— 

Even without Coulomb interaction U this model is not 
exactly solvable in general, thus approximations have to 
be made. For U = 0, a first approach can consist of 
a perturbative treatment of the exchange interaction J, 
leading to the well-known effective Rudcrmann-Kittel- 
Kasuya-Yoshida (RKKY) interactio n 23 ' 24 ^ 5 with a char- 
acteristic, dimension-dependent dependence on distance. 
Although generally accepted as an at least proper ansatz 
for a qualitative discussion of magnetic properties of 
models like the KLM, it has not yet been studied in detail 
how well this approximation works for increasing J. Fur- 
thermore, for the antiferromagnetic Kondo lattice model, 
where heavy Fermion physics can play an essential role, 
or in the presence of additional correlations in the band 
states, the validity of the use of the RKKY arguments is 
far from clear. Thus one aspect of the present paper is to 
investigate to what extent the RKKY exchange indeed 
leads to a reasonable description of the low-temperature 
properties of the KLM. 

As approximation to study the KLM while leav- 
ing as much of the local correlations induced by both 
the Coulomb interaction U and the exchange J intact 
as possible we use the Dynamical Mean Field Theory 
(DMFT))2£ mapping the lattice model onto an effective 
impurity problem. In former treatments of the model |T]), 



especially within DMFT ) 27 i 28 i 29 classical spins were as- 
sumed to avoid the sign problem of the quantum Monte- 
Carlo treatment of the effective impurity problem em- 
ployed there. The importance of a fully quantum me- 
chanical treatment of the local spins even in impurity 
problems has been addressed by Peters and Pruschke^ 
and the effects of quantum spins in the KLM by Kienert 
and Nolting— . To achieve such a fully quantum mechan- 
ical treatment we use Wilsons's Numerical Rcnormaliza- 
tion Group (NRG)^1^ as impurity solver. The NRG 
can handle the whole interaction and temperature regime 
and is also able to calculate Green's function even in or- 
dered phases. We use a recent improved technique for 
calculating Green's function a 35 ' 36 i 37 within NRG. There- 
fore we are able to treat a large bandwidth of values of 
both on-site interactions J and U. 

The paper is organized as follows. In the next section 
we discuss the form of the perturbative RKKY exchange 
in the limit spatial dimension D — ► oo appropriate for 
the DMFT. In section III we present our results for the 
magnetic phase diagrams of the model (JTJ) at half-filling 
and finite doping and different values of U and discuss 
their dependence on the sign and magnitude of J. A 
summary will conclude the paper. 



II. RKKY-INTERACTION IN THE LIMIT 

D — > oo 

Conventionally the interaction between localized 
spins in metals is analyzed in terms of the RKKY 
interaction i 23 i 24 i 25 This effective exchange interaction 
shows a characteristic dependence on distance R and lo- 
cal exchange coupling J of the form 



2 cos(k F R) 

Jrkky k j -jk^Rr 



(2) 



where for is the Fermi momentum of the host and a some 
positive, dimension-dependent number. This distance- 
dependence is due to the sharp Fermi surface in metals 
and, strictly speaking, valid only in the limit J — > 0. 
Note that, because Jrkky oc J 2 , the sign of J does 
not matter. Furthermore, it is a priori not evident, how 
the distance-dependence looks like in the limit spatial 
dimension D — > oo, which sets the framework for the 
construction of the DMFT. In particular, in this limit 
there is no proper definition of for, which controls the 
dependence of Jrkky on the occupancy of the band. 
Thus, in order to be able to study to what extent the 
RKKY interaction controls the magnetic properties of 
the KLM in DMFT, we will derive this interaction and 
especially its dependency on the filling for the limit D — > 
oo in the following. 

To this end we consider the extended Hubbard model 
(P) on a hypercubic lattice. In the paramagnetic phase 
with (Si) = (si) = 0, an effective Hamiltonian for the 
local spins can be calculated pcrturbativcly in J by for- 
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mally tracing out the band states. The resulting expres- 
sion in lowest order in J reads 

3 

h cS » - ^ 2 E E / dT s " s " ( s ? m s j t°] ) j=o 
= -^EE^w > 

ij ct— 1 

where xTf denotes the static susceptibility of the bare 
Hubbard model. In the paramagnetic case, xff does not 
depend on a, and we are free to evaluate it for a = 3, for 
example. This leads to an effective spin-spin interaction 

Jrkky {Ri - Rj) = J Xij 



= ■^£ e *' ( *"* ) x"(9) (3) 

1 

Let us evaluate the latter sum for a model with nearest- 
neighbor hopping in the limit D — ► oo. As has been 
shown by Miiller-Hartmann^ the ^-dependence then en- 
ters only via 



v(q) ■= ^E cos ( fta ) 

i=i 

where a is the lattice parameter. 
Due to inversion symmetry, we have 



l 

±_Y^e^ R >- R ^x zz (q) = 2^E C0S (« ' - R i))x zz (q) = J dx Qij (x) X zz (x) 
i i _i 

I 



with 

gij (x) := 2l [cos (q ■ (Ri - Rj)) 5(x - r,(q))] . 
1 

For a nontrivial result in the limit D — > 00 we need D ■ 
Jrkky to be finite. We thus will evaluate 

1 

D / dx QlJ {x)x zz {x) 



directly. Following again the arguments by Miiller- 
Hartmann, we first rewrite 



6{x- V (q)) = J 

— OO 

and obtain 



7T OO 

D Qij (x)=2D J J ^-v(i))s cosiq . iRi _ Rj)) 

— 7T —OO 



Expanding the exponential in terms of 77(g) and observing for i =/= j 
that 



■cos (g- (iZj - -Rj)) 



— cos^-i^gcos^-a) = ^_ flj | >a , -iS ]Ri . RjUa J ^s + oU 

— OO 

-5 — 5(x) O (—} 

and so on for terms involving rj(q) m for m > 1, we obtain |-Ri-fij|,o ^ I J \DJ 
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With this result we find in the limit 
RKKY exchange 



D 



oo for the 



D ■ Jrkky (Ri — Rj ) — J 



(x) 



(4) 



x=0 



Note that the right hand side of ((4]) is already correctly 
scaled to obtain nontrivial results in the limit D — > oo. 
Furthermore, the RKKY exchange acts only on nearest 
neighbors, the oscillatory structures arising from Fermi 
surface singularities in finite dimensions are absent in 
D — ► oo. Nevertheless, the dominant nearest-neighbor 
exchange constant which through its sign determins the 
type of order can be obtained. 

For U = 0, the susceptibility is given by the simple 
bubble, which can be evaluated exactly albeit only nu- 
merically, leading to the behavior for Jrkky depicted 
in Fig. [TJ Note that the sign change from J < (afm) 
to J > (fm) appears for n < 0.5. We will discuss this 
figure in connection with numerical results in the next 
sections. 
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1: -j^X zz ( x )\x=o for U = as function of the occupancy 



For finite U > 0, a similar evaluation is not possible, 
because the susceptibility entering now is the full lattice 
susceptibility. While Xzz(x = 0) is the local susceptibil- 
ity, which can be obtained from the effective impurity 
problem directly, its derivative determining Jrkky in- 
volves neighboring x values, which cannot be calculated 
within NRG. We content ourselves here by noting that fi- 
nite Coulomb correlations typically lead to a more asym- 
metric distribution of spectral weight, which is known to 
tend to enhance ferromagnetic correlations. We thus ex- 
pect that the root of Jrkky will shift to larger values of 
the occupancy n for finite U . 



III. RESULTS 

In this section wc present our results for the magnetic 
phases of the extended KLM |T]). Up to now the lo- 



cal spin S was completely arbitrary and can in prin- 
ciple take any value. Although it is surely interesting 
to study the effect of increasing spin quantum number 
on the results^ we restrict the discussion to the case 
5 = 1/2 here. 

The calculations were done for a Bcthe lattice em- 
ploying a bipartite subdivision to accommodate antifer- 
romagnctic order. 26 ! 39 We used a Bethe lattice instead of 
a hypcrcubic one for computational reasons. However, 
we did not find any significant differences between the 
two lattices in test calculations. As discretization pa- 
rameter for the NRG we used A = 2 and typically kept 
1000 . . . 2000 states in each NRG step. As our unit of 
energy we choose the bandwidth W = At. 



A. Antiferromagnetism at half filling 

Let us begin by examining the magnetic order at half- 
filling, n = 1. Clearly Fig. [T] states that the interac- 
tion between the localized spins is negative which is sup- 
posed to result in antiferromagnetic order. In Fig. [2] are 
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FIG. 2: (color online) Antiferromagnetic polarization at half- 
filling and T — as function of the coupling J between spins 
and band electrons. Circles were calculated for U/W = 1 /200, 
corresponding approximately to U = 0, squares for U — W . 

collected results for the polarization (s z ) = — njj 

of the band electrons and (S z ) of the local spin. One 
can see in Fig. [2] two curves corresponding to U f=a 
(circles) and U = W (squares). For large negative J 
there is no magnetic order. The system forms a Kondo- 
insulator, locally quenching all moments. Within DMFT 
we find for the critical value of the coupling at U = 
J c /t w — 0.3W/t = —1.2 which is consistent with long- 
known results4£ With increasing U, this value is shifted 
to somewhat smaller absolute values \J\. 

For J > J c antiferromagnetic order can emerge before 
all moments have been quenched locally. Quite gener- 
ally, we observe that the polarization (s z ) of the band 
states and (S z ) of the local spins are opposite in sign 
for J < 0. Around the critical value J c the resulting 
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total polarization (s z + S z ) thus is very small, although 
both contributions can already have rather sizable values. 
This result is quite interesting, in particular in view of 
the so-called small-moment antiferromagnetism observed 
in several rare-earth compounds. 

Obviously, for U — there is no antifcrromagnctic 
order at J = 0; however, even for very small |J/W| -C 1 
the local spins are almost fully polarized, (S z ) « 0.5. 
On the other hand, the polarization of the conduction 
electrons (s z ) goes smoothly through zero. We can thus 
identify this range of J values as corresponding to the 
RKKY regime, where the local spins are fully polarized 
and the band electrons show a polarization proportional 
to the "effective field" ~ J(S Z ) provided by the local 
spins. In this region we also observe that the magnetic 
properties of the system are roughly independent of the 
sign of J, as predicted by RKKY. 

On the other hand, the behavior in the vicinity of J c 
cannot be understood in terms of RKKY any more, al- 
though | J c \ = 0.3W is still significantly smaller than W 
and one might expect the perturbational arguments to be 
still valid. However, the physical properties are radically 
different for J < and J > 0, as no critical point exists 
for J > and the model always is in the fully polarized 
state. 

For U > there is antifcrromagnetic order even at 
J = 0, which represents the pure Hubbard modeL 26 i 39 
As mentioned before, for J < spin and electron polar- 
izations are antiparallcl while for J > both orientations 
arc parallel. Due to to numerical errors we were not able 
to determine the behavior as J — > 0. Even with as many 
as 4000 states kept in each NRG step DMFT+NRG cal- 
culations did not converge but showed strong fluctuation 
for | J/W\ < 0.04. The results obtained nevertheless sug- 
gest a jump at J = 0. Again, as J — > 0~, the net polar- 
ization (s z + S z ) is strongly reduced from the almost full 
values for each individual part. 




Fig. [3] shows the majority spectral functions for 
U/W = 1/200 and different J. All calculations show 
a gap at w = 0, characteristic for an insulating behav- 
ior. The solid line represents the results for a large 
antifcrromagnetic coupling, i.e. in the Kondo insulat- 
ing state. The spectral function agrees with the typi- 
cal result for paramagnetic phase of a Kondo-insulator. 
The remaining lines correspond to J > J c . The dotted 
and dashed lines correspond to the same value \J\, how- 
ever with opposite signs. The general shape is that of 
a weak-coupling antiferromagnet^ with the characteris- 
tic square-root van-Hove singularities at the gap edges. 
Note that for J < we find a shallow shoulder for uj > 
which quite likely is the remnant of a quasiparticle peak 
due to the Kondo effect expected in this regime^ This 
is again a sign, that already for moderate values of J the 
physical properties are significantly different for J < 
and J > 0, i.e. cannot be understood by pure RKKY 
physics. 

The dotted-dashed line in Fig. [^corresponds to a large 
ferromagnetic coupling. Here, life-time effects due to a 
large self-energy already substantially dominate the spec- 
trum. 



B. Doped model at U — 
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FIG. 3: (color online) Spectral function for the majority spin 
at half filling for U/W = 1/200 and temperature T/W = 
2.3 ■ 10~ 4 . 
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FIG. 4: (color online) Magnetic phase diagrams for (7 = 
and | J\ = W/2 as function of filling and temperature. The 
upper (lower) panel shows the results for ferromagnetic (anti- 
ferromagnetic) coupling. Lines are meant as guide to the eyes. 
For ferromagnetic coupling a magnetic phase exists between 
0.45 < n < 0.7, whose nature could not be determined (see 
text). The number of states kept during NRG calculation was 
increased at phase boundaries. 

Figure 5] shows the magnetic phases of the Kondo- 
lattice model at U — as function of filling n and tem- 
perature calculated within DMFT for a Kondo coupling 
\J\ = W/2. Note that this value is already larger than 
the critical value J c at half filling, thus we expect that 
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for J < the physics is governed by local quenching of 
the moments close to half filling. 

The upper panel in Fig. [4] shows the model with fer- 
romagnetic coupling, while the lower one displays the 
results for antifcrromagnctic coupling. In the ferromag- 
netic Kondo-lattice model we find an antiferromagnetic 
phase around half-filling n > 0.7, and a ferromagnetic 
phase for n < 0.45. The location of phase boundary 
of the ferromagnet at T = appears to agree very well 
with RKKY prediction from Fig. [T] where the coupling 
between the localized spins within DMFT becomes fer- 
romagnetic for n < 0.45. We have focused in our work 
on homogeneous ferromagnetic-, antiferromagnetic Neel- 
and paramagnetic states. Between 0.45 < n < 0.7 
calculations showed also significant magnetic response. 
But neither a ferromagnetic nor an antifcrromagnctic 
Neelstate could be stabilized. The real nature of this 
phase could not yet be clarified in our calculations, but 
one might argue that one should expect some type of 
incommensurate order. In fact, the magnetic phase dia- 
gram of the Kondo-lattice model has been analyzed by a 
number of other authors, 31 i 40 ' 42 i 4 V 4 i 45 ' 46 i 47 i 4i M 9 In addi- 
tion to the conventional homogeneous ferromagnet and 
Neel antifcrromagnct other magnetic phases like incom- 
mensurate, chiral and short range ordered phases were 
analyzed for 0.5 < n < 1.0 in these studies and found 
to be stable in the regime where our calculations fail to 
converge. 
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FIG. 5: (color online) Occupancy (upper panel) and polariza- 
tion (lower panel) for the antiferromagnetic solution around 
half filling as function of chemical potential /j, for different 
temperatures. Other parameters as in Fig. [4] 

For the antifcrromagnet close to n = 1, several 
author o 31 i 43 ' 44 i 45 ' 48 also reported phase separation be- 
tween antiferromagnetism at half-filling and phases away 
from half-filling. Figure [5] shows our results for the oc- 
cupancy and polarization in the antiferromagnetic phase 
as function of the chemical potential [i for different tem- 
peratures. We see that with decreasing temperature the 
slope of the curves for n — > 1.0 get larger and larger. 



However, even for the lowest temperature we could sta- 
bilize every occupation number n < 1.0, even if there is 
a large slope for n — > 1.0. Increasing the numerical accu- 
racy by for example keeping more states per NRG step 
also tend to increase the slope further. However, we did 
never observe a clear jump of the occupation at a criti- 
cal value of the chemical potential. Note, however, that 
we cannot rule out phase separation at T = from our 
numerical results. 

The lower panel in figure 0] shows the phase diagram 
of the antiferromagnetic Kondo-lattice model. Here our 
results yield only evidence for a ferromagnetic phase. 
In contrast to the ferromagnetic Kondo-lattice model, 
this phase reaches to significantly higher occupations and 
temperatures; the antifcrromagnctic coupling seems to 
stabilize ferromagnetic order. This behavior as well as 
the rather large critical value for n at T = cannot be 
explained within the RKKY picture. Furthermore, the 
quite likely incommensurate phase neighboring the ferro- 
magnet for J > is absent here. For n — ► we could not 
stablilize a magnetic phase any more due to numerical 
problems connected with the small filling. 
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FIG. 6: (color online) Polarization of the ferromagnetic solu- 
tions. The upper (lower) panel corresponds to ferromagnetic 
(antiferromagnetic) coupling. Other parameters as in Fig. [4] 

The previous observations are further substantiated by 
the results in Fig. [51 where we show the ferromagnetic po- 
larizations for ferromagnetic (upper panel) and antiferro- 
magnetic (lower panel) coupling parameter. The circles 
correspond to a very low temperature, squares to a tem- 
perature near the transition into the paramagnetic state. 
Full lines represent (s z ), dashed lines (S z ). For J > and 
temperatures above the critical temperature for the "in- 
commensurate" phase, both (s 2 ) and (S z ) show a similar 
behavior as function of doping. In particular, they van- 
ish smoothly as n — ► n c , indicating a second order phase 
transition. For very low T, on the other hand, (S z ) re- 
mains almost constant at the fully polarized value until 
we reach the phase boundary to the "incommensurate" 
phase at n c ps 0.45. At the same time, (s z ) monotoni- 
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cally increases with filling up to n c . Although we cannot 
determine the true structure of the phase for n> n c , its 
existence seems to be connected to the peculiar behav- 
ior of {s z } rather than with (S z }, i.e. is primarily driven 
by the itinerant electrons. We are currently further in- 
vestigating this phase transition by including magnetic 
structures other than the Neel state in our calculations. 

For antifcrromagnctic coupling (lower panel of Fig. 
[6]) we find a critical n c increasing with temperature, 
where the ferromagnetic metal becomes a Kondo insu- 
lator. Note that for T/W = 4 • 10~ 5 we are effectively in 
the ground state. Further lowering the temperature does 
not change the phase diagram any more. Both (s z ) and 
(S z ) appear to vanish smoothly, too, thus again indicat- 
ing a second order phase transition. 

Choosing \J\ < \J C \ does not modify the general struc- 
ture of the phase diagram for ferromagnetic coupling. 
The results for antiferromagnetic coupling, of course, do 
change, in particular close to half filling, where now an 
antiferromagnetic phase appears, too. While the overall 
structure of the phase diagram is now similar to the case 
J > in Fig. |4] upper panel, the values for are 
still enhanced with respect to the RKKY prediction. 



C. Finite Coulomb interaction U = W/2 
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FIG. 7: (color online) Magnetic phase diagrams for U = W/2 
and | J = W/2 as function of filling and temperature. The 
upper (lower) panel shows results for ferromagnetic (antifer- 
romagnetic) coupling. Lines are meant as guides to the eye. 
For ferromagnetic coupling a magnetic phase exists between 
0.5 < n < 0.8, whose nature could not be determined (see 
text). The number of states kept during NRG calculation 
was increased at phase boundaries. 

Let us now turn to the model with finite on-site 
Coulomb interaction. We present results for fixed U = 
\J\ = W/2, changing these parameters modifies the 
observations quantitatively but not qualitatively. The 
phase diagrams for ferromagnetic (upper panel) and an- 



tifcromagnctic (lower panel) coupling are shown in Fig. 
[7] As for [7 = 0, we find both an antifcrromagnctic 
phase around half filling and a ferromagnetic phase for 
larger doping in the ferromagnetic Kondo-latticc model. 
The critical values have changed to ric^ F ^ « 0.8 and 
n\f M ^ w 0.53 for T = 0, indicating that local correlations 
additionally stabilize ferromagnetic order but destabilize 
antifcrromagnctic for J > 0, in accordance with the an- 
ticipated variation of Jrkky with U. 

For the ferromagnetic Kondo lattice model we again 
find between the ferromagnetic phase for n < 0.53 and 
the antiferromagnetic for 0.8 < n < 1 a magnetic phase 
whose nature could not be determined for the same rea- 
sons as for U — 0. Guided by our observations for U = 
and a comparison to other studies by other groups, one 
can again assume the appearance of an incommensu- 
rate magnetic phase like a spin wave or a chiral phase. 
From quantum Monte-Carlo studies^ for the conven- 
tional Hubbard model an incommensurate phase can in- 
deed be anticipated in this parameter regime. 
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FIG. 8: (color online) Occupancy (upper panel) and polariza- 
tion (lower panel) for the antiferromagnetic solution around 
half filling as function of chemical potential /i for different 
temperatures, chemical potential fi. Other parameters as in 
Fig. [7] 

Figure[8]shows occupation and polarization as function 
of the chemical potential [i for different temperatures. 
For high-temperatures, squares and triangles, there is a 
smooth crossover from half-filling to lower occupancies, 
as for U = (cf. Fig. However, in contrast to the 
indecisive results at U = 0, we observe clear evidence for 
phase separation at T/W = 4 ■ 10 -5 , circles, this time. 
There exists a critical value fJ-critj where both occupation 
and polarization jump from their values at half-filling to 
a smaller occupancy and polarization. We therefore have 
phase separation between the antiferromagnetic insulator 
at half-filling and an antiferromagnetic metal with a lower 
occupancy. Note that one also finds phase separation 
in the conventional Hubbard model, which here however 
occurs between an antifcrromagnctic insulator with n = 1 
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and a paramagnet with n < n 
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FIG. 9: (color online) Polarization of the ferromagnetic solu- 
tions. The upper (lower) panel corresponds to ferromagnetic 
(antiferromagnetic) coupling. Other parameters as in Fig. [Jj 
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IV. SUMMARY 



The dependence of the polarization on n and T in the 
ferromagnetic phase for the present values of U and \J\ is 
shown in Fig. [5] for J > (upper panel) and J < (lower 
panel). As in Fig. [6l full lines represent (s z ) and dashed 
lines (S z ). For high temperatures (squares) we again find 
second order phase transitions, indicated by the smooth 
simultaneous vanishing of both (s z ) and (S z ), while for 
low T (circles) and J > the peculiar behavior at the 
phase boundary to the "incommensurate" phase already 
seen in Fig. [5] is observed. 

A markable difference to the behavior for U = col- 
lected in Fig. [6] can be observed for antiferromagnetic 
coupling J < 0. Here, the polarization now vanishes dis- 
continously at n^f for T — > 0, indicating a first order 
phase transition in contrast to the second order type ob- 
served for (7 = 0. Thus, while the results at U = would 
predict a quantum critical point at ni FM \ such a sce- 
nario is obvioulsy destroyed by Coulomb correlations in 
the conduction band, at least within DMFT. 

Note that for the bare Hubbard model no ferromag- 
netic phase exists in this parameter regime for a bipar- 
tite lattice. It is thus clearly the spin-electron interaction 
which leads to ferromagnetism. Our results are consistent 
with earlier studies by other groups, who however treated 
the electron-electron interaction in an effective medium 
approach , 51 i 52 ' 53 i 54 while we are able to treat the on-site 
correlations exactly within DMFT and NRG. Especially 
our Curie temperatures are located in the same range as 
in these earlier studies. 

A possible reason for the tendency of both finite U 
and negative J to stabilize ferromagnetism is the strong 
asymmetry induced by U and J 55 i 56 . As can be seen 
from Fig. [TUl this effect is most pronounced for J < 0. 
For J > 0, a weaker, but nonetheless significant redistri- 
bution of spectral weight is observed, which can explain 
the increase of n^' M ^ here, too. 



In this paper we discussed the magnetic properties of 
the extended Kondo lattice model, where in addition to 
an exchange coupling between itinerant states and lo- 
calized spins the band electrons are subject to a local 
Coulomb interaction. Such a model is expected to de- 
scribe the properties of e.g. transition metal compounds, 
where part of the d-states are localized in the crystalline 
environment due Coulomb correlations and form a local 
spin which couples to the remaining, usually itinerant 
d-states either due to Hund's exchange or hybridization. 

These materials typically have a rich magnetic phase 
diagram, and it is an interesting question, how these 
two interactions cooperate or compete in driving differ- 
ent magnetic phases. A particularly interesting question 
in this respect is, to what extent, qualitatively and quan- 
titatively, the concept of the RKKY approximation of an 
effective interaction between the local spins mediated by 
the itinerant conduction states holds. 

We treat the model within DMFT, using NRG as solver 
for the effective impurity problem. The calculations were 
done for a bipartite lattice, allowing for homogeneous fer- 
romagnetism and Neel antiferromagnetism. At half fill- 
ing, we find an antiferromagnetic phase, as expected, for 
all ferromagnetic Kondo couplings J > 0. For antiferro- 
magnetic coupling J < 0, there exists a critical J c , where 
the system becomes a paramagnetic Kondo insulator. 1 ® 

For J > J c , this antiferromagnetic phase prevails for 
finite doping up to a critical filling ric depending on 
J and the local Coulomb interaction U . For finite U and 
J > 0, we in addition find phase separation between the 
antiferromagnetic insulator with n = 1 and an antifer- 
romagnetic metal with n < 1. For U = we find no 
convincing evidence that such a phase separation exists, 
too, but the numerical results are not sufficiently clear to 
rule it out cither. For J < J c no antiferromagnetic phase 
exists at all. 
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Below a second critical filling rv c ' a homogeneous 
ferromagnet is found. Quite interestingly, the extent and 
stability of this phase is increased by both a finite U and 
an antiferromagnetic J. The reason for this stabiliza- 
tion can quite likely be traced back to the introduction 
of a strong asymmetric redistribution of spectral weight 
in the electronic spectral function. An interesting aspect 
in connection with quantum critical points in e.g. ferro- 
magnetic heavy fermion compounds is that for antiferro- 
magnetic Kondo exchange a finite Coulomb correlation 
among the conduction states has the tendency to change 
the order of the phase transition at ni FM ' from second 
order without correlation to first order with correlations. 

We find that for ferromagnetic Kondo coupling J the 
predictions by RKKY can be used to at least qualita- 
tively account for the different phases. For moderate and 
large J < 0, however, RKKY cannot even qualitatively 



predict the phases correctly, underestimating the ferro- 
magnetic phase grossly. Moreover, the variations of the 
polarizations of the local spin and bandelectrons seem 
to follow RKKY for very small |J|, but already deviate 
substantially for moderate values. 
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